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Abstract 

We report results of quark masses in quenched lattice QCD with the Kogut- 
Susskind fermion action, employing the Reguralization Independent scheme 
(RI) of Martinelli et al. to non-perturbatively evaluate the renormalization 
factor relating the bare quark mass on the lattice to that in the continuum. 
Calculations are carried out at (3 = 6.0, 6.2, and 6.4, from which we find 
m^ s (2GeV) = 4.23(29)MeV for the average up and down quark mass and, 
with the d> meson mass as input, m^ s (2GeV) = 129(12)MeV for the strange 
mass in the continuum limit. These values are about 20% larger than those 
obtained with the one-loop perturbative renormalization factor. 

PACS number(s): 11.15.Ha, 12.38.Gc, 12.15.Ff 
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The values of quark masses are fundamental parameters of the Standard Model which 
are not directly accessible through experimental measurements. Lattice QCD allows their 
determination through a calculation of the functional relation between quark masses and 
hadron masses. For this reason a number of lattice QCD calculations have been carried out 
to evaluate quark masses, employing the Wilson, clover or Kogut-Susskind (KS) fermion 
action 

An important ingredient in these calculations is the renormalization factor relating the 
bare lattice quark mass to that in the continuum. While perturbation theory is often used 
to evaluate this factor, uncertainties due to higher order terms are quite significant in the 
range of the QCD coupling constant accessible in today's numerical simulations. A non- 
perturbative determination of the renormalization factor is therefore necessary for a reliable 
calculation of quark masses, and effort in this direction has recently been pursued for the 
Wilson and clover fermion actions |@-|5|. 

The need for a non-perturbative determination of the renormalization factor is even 
more urgent for the KS action since the one-loop correction ||[7| is as large as 50% in 
present simulations. In this article we report a study to meet this need |8],|9| : we calculate 
the renormalization factor of bi-linear quark operators for the KS action non-perturbatively 
using the Reguralization Independent scheme (RI) of Ref. 0] developed for the Wilson/clover 
actions. The results for the scalar operator, combined with our previous calculation of bare 



quark masses ||10|| , lead to a non-perturbative determination of the quark masses in the 
continuum limit. 

In the RI sheme, the renormalization factor of a bi-linear operator O is obtained from 
the amputated Green function, 

r (p) = sipYHWip) o fomsip)- 1 (l) 

where the quark two-point function is defined by S{p) = (O|0(p)0(p)|O). The quark field <fi(p) 
with momentum p is related to the one-component KS field x{ x ) by 4>a (p) = S y ex P(~ W • 
y)x(v + aA), where y^ = 2an M , p^ = lim^jiaV) with — L/4 < < L/4 and = 0, 1. 
Bi-linear operators have a form 

o= E 4> a A{y)iis®&) AB u a A b B (y)<P b B (y) (2) 

yABab 



where (75 <8> £f) refers to Dirac (75) and KS flavor structure 0, and the indices a and 
b refer to color. The factor U^ B {y) is the product of gauge link variables along a minimum 
path from y + aA to y + aB. We note that Uab{3J) is absent for scalar and pseudo scalar 
operators as these operators are local. 

The renormalization condition imposed on Tq(p) is given by 

Z^\p)-Z^p)=Tt[pIt (p)] (3) 



where Pq = (75 <S> is the projector onto the tree-level amputated Green function. The 
wave function renormalization factor Z^(p) can be calculated by the condition Zy(p) = 1 
for the conserved vector current corresponding to (7^ <8> /). Since the RI scheme explicitly 
uses the quarks in external states, gauge fixing is necessary. We employ the Landau gauge 
throughout the present work. 
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The relation between the bare operator on the lattice and the renormalized operator in 
the continuum takes the form, 



^ms(m) = U m (^p)Z^(p)/Z$(p)0 (4) 

where U^(/i, p) is the renormalization-group running factor in the continuum from momen- 
tum scale p to fi. We adopt the naive dimensional regularization (NDR) with the modified 
minimum subtraction scheme (MS) in the continuum. The factor Z^(p) provides matching 
from the RI scheme to the MS scheme. These two factors are calculated perturbatively in 
the continuum. For our calculation of the quark mass we apply the relation @ in the scalar 
channel in the chiral limit, i.e, 1/Z^ = 

Our calculations are carried out in quenched QCD. Gauge configurations are generated 
with the standard plaquette action at (3 = 6.0, 6.2, and 6.4 on an 32 4 lattice. For each /3 we 
choose three bare quark masses tabulated in Table | where the inverse lattice spacing 1/a 
is taken from our previous work [[Kj. 

We calculate Green function for 15 momentum in the range 0.038553 < (ap) 2 < 1.9277. 
Quark propagators are evaluated with a source in momentum eigenstate. We find that the 
use of such a source results in very small statistical errors of 0(0.1%) in Green functions. 

The RI method completely avoids the use of lattice perturbation theory. We do not have 



to introduce any ambiguous scale, such as q* [11|, to improve on one- loop results. An impor 



tant practical issue, however, is whether the renormalization factor can be extracted from 
a momentum range Aqcd <C p <C 0(1/ a) keeping under control the higher order effects in 
continuum perturbation theory, non-perturbative hadronization effects, and the discretiza- 
tion error on the lattice. These effects appear as p dependence of the renormalization factor 
in (|j), which should be absent if these effects are negligible. 

In Fig. H we compare the scalar renormalization factor Zf l (p) with that for pseudo scalar 
Zp ll (p) for three values of bare quark mass am at (3 = 6.0. From chiral symmetry of the 
KS fermion action, we naively expect a relation Z^ l (p) = Zp l (p) for all momenta p in the 
chiral limit. Clearly this does not hold with our result toward small momenta, where Zp l (p) 
rapidly increases as m — > 0, while Z^ l (p) does not show such a trend. 

To understand this result, we note that chiral symmetry of KS fermion leads to the 
following identities between the amputated Green function of the scalar T s (p), pseudo scalar 
Tp(p), and the quark two-point function S(p) _1 : 

r*(p) = ^(pr 1 ( 5 ) 



Tp(p) = ^[(js®tis)S(p)- 1 + S(p)- 1 (K®&)] (6) 
We also find numerically that the quark two-point function can be well represented by 

S(p)- 1 ~ 53 (^¥7M ( p )zC» + M(p) (7) 



with two real functions C^(p) and M(p), where £t(p) = cos(ap fl ) — ^(7^75 <8> £^£5) sin(ap M ). 
From d6|), (0), and @ we obtain the relations between the renormalization factors and M(p), 

Zf(p)-Z lP (p)=dM(p)/dm (8) 
Zf(p)-Z^(p)=M(p)/m (9) 
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In Fig. [J M(p) in the chiral limit obtained by a linear extrapolation in m is plotted. It 
rapidly dumps for large momenta, but largely increases toward small momenta. Combined 
with (|9]) this implies that Zp l (p) diverges in the chiral limit for small momenta, which is 
consistent with the result in Fig. |l|. 

The function M(p) is related to the chiral condensate as follows : 

(10) 

A non-vanishing value of M(p) for small momenta would lead to a non-zero value of the 
condensate. Therefore the divergence of Zp l (p) near the chiral limit is a manifestation of 
spontaneous breakdown of chiral symmetry; it is a non-perturbative hadronization effect 
arising from the presence of massless Nambu-Goldstone boson in the pseudo scalar channel. 

While we do not expect the pseudo scalar meson to affect the scalar renormalization 
factor Zg- l (p), as indeed observed in the small quark mass dependence seen in Fig. U, the 
above result raises a warning that Z<p(p) may still be contaminated by hadronization effects 
for small momenta. 

In Fig. |3] we show the momentum dependence of Z m (fj,,p, 1/a) = Uy^(n,p)Z^i(p)Z^ l (p) 
which is the renormalization factor from the bare quark mass on the lattice to the renormal- 
ized quark mass at scale // in the continuum. Here we set \i = 2GeV and use the three-loop 
formula |L2| for Uy^ and Z^j. While Z m (fi,p,l/a) should be independent of the quark 



momentum p, our results show a sizable momentum dependence which is almost linear in 
(ap) 2 for large momenta (filled symbols in Fig. [|). 

For small momenta we consider that the momentum dependence arises from non- 
perturbative hadronization effects on the lattice and the higher order effects in continuum 
perturbation theory. It is very difficult to remove these effects from our results. 

Toward large momenta, however, these effects are expected to disappear. The linear 
dependence on (ap) 2 , which still remains, should arise from the discretization error on the 
lattice, i.e., 

Z m (fi,p, 1/a) = m m (fi)/m + (ap) 2 Z n + 0(a 4 ) (11) 

with the constant Z-r corresponding to the mixing to dimension 5 operators on the lattice. 

This relation implies that, if we take a continuum extrapolation of Z m ([i,p, l/a)m at a 
fixed physical momentum p, the discretization error in Z m is removed. This procedure also 
removes the a 2 discretization error in the lattice bare quark mass m itself reflecting that in 
hadron masses. 

The momentum p should be chosen in the region where the linear dependence on (ap) 2 
is confirmed in our results. This region starts from a similar value of p 2 ~ 3GeV 2 for the 
three (3 values, and extends to p 2 ~ 1.9/ a 2 , the highest momentum measured. Hence we are 
able to use only a rather narrow range 3GeV 2 < p 2 < 6.6GeV 2 , the upper bound dictated 
by the value of 1.9/a 2 for the largest lattice spacing at /3 = 6.0. 

In Fig. f| we show the continuum extrapolation for the averaged up and down quark 
mass at fi = 2GeV. Filled circles are obtained for p = 1.8GeV and squares for p = 2.6GeV 
for which the value of Z m is obtained by a linear fit in (ap) 2 employing the filled points in 
Fig. |||. The bare quark mass |T(| is determined by a linear extrapolation of pseudo scalar 
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meson mass squared in the Nambu-Goldstone channel (75 ® £5) and that of vector meson 
mass in the VT channel (7$ ® £j) to the physical point of tt and p meson masses. 

We observe that the continuum extrapolation completely removes the momentum depen- 
dence of the quark mass at finite lattice spacings. Furthermore the values are substantially 
larger than those obtained with one-loop perturbation theory (open circles for q* = 1/a and 
squares for q* = it /a ). Making a linear extrapolation in a 2 , our final result in the continuum 
limit is 



ni 



MS 
ud 



(2GeV) = 4.23(29)MeV. (12) 



where we adopt the value for p = 2.6GeV since this is the largest momentum accessible 
and the momentum dependence is negligible. This value is about 20% larger than the 
perturbative estimates : 3.46(23)MeV for q* = 1/a and 3.36(22)MeV for q* = ir/a. We 
collect the values of renormalization factor and quark masses in Table and [TT1|. 
Applying our renormalization factor to the strange quark mass, we obtain 



m 



A1S 



(2GeV) = 106.0(7.1)MeV for m K (13) 
= 129(12)MeV for (14) 



where we use K or meson mass to determine the bare strange mass. Results from pertur- 
bative estimation are given in Table |I^. 

The CP-PACS Collaboration recently reported the results fT3J m^(2GeV) = 4.6(2)MeV, 



m 



MS 



(2GeV) = 115(2)MeV(mif) and 143(6)MeV(m^) from a large-scale precision simulation 
of hadron masses with the Wilson action. Our values are 10% smaller, which may be due 
to the use of one-loop perturbative renormalization factor in the CP-PACS analysis. 
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Accelerator Research Organization (KEK), and also in part by the Grants- in- Aid of the 
Ministry of Education (Nos. 08640404, 09304029, 10640246, 10640248, 10740107, 10740125). 
S.K. and S.T. are supported by the JSPS Research Fellowship. 
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FIG. 1. The renormalization factor for the scalar Zg J (p) and the pseudo scalar Zp I (p). 
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FIG. 2. M(p) in the chiral limit. 
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FIG. 3. The ratio Z m (fi,p, 1/a) at jj, = 2GeV. For each [3 the filled data points are used for 
linear interpolation in (ap) 2 . 
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FIG. 4. The final results of the light quark mass m^ s (2GeV) . 



10 



TABLES 



(3 


l/o(GeV) 




am 




# conf. 


6.0 


1.855(38) 


0.010 


0.020 


0.030 


30 


6.2 


2.613(87) 


0.008 


0.015 


0.023 


30 


6.4 


3.425(72) 


0.005 


0.010 


0.020 


50 



TABLE I. Run parameters. 



NP. P. 





p = 1.8GeV 


p = 2.6GeV 


q* = 1/a 


q* = ir/a 


6.0 


2.735(16) 


2.509(18) 


1.867 


1.776 


6.2 


2.5451(91) 


2.449(10) 


1.877 


1.800 


6.4 


2.385(12) 


2.339(12) 


1.871 


1.804 



TABLE II. The non-perturbative (NP) and the perturbative (P) renormalization factors of 
the quark mass Z m (/j,,p, 1/a) at ji = 2GeV. 



& 


am ud 

(io- 4 ) 


NP.(MeV) 
p = 1.8GeV p = 2.6GeV 


P.(MeV) 
g* = 1/a q* = ir/a 


6.0 


9.72(40) 


4.93(23) 


4.52(21) 


3.37(15) 


3.20(15) 


6.2 


6.86(45) 


4.56(34) 


4.39(32) 


3.36(25) 


3.23(24) 


6.4 


5.38(23) 


4.39(21) 


4.31(21) 


3.45(16) 


3.32(16) 


a 2 ^0 




4.17(30) 


4.23(29) 


3.46(23) 


3.36(22) 



TABLE III. Final results for m£j(2GeV) obtained with the non-perturbative (NP) and the 
perturbative (P) renormalization factors. The lattice bare quark mass am U( i is also listed. 



p 


am s 


NP.(MeV) 


P.(MeV) 




(io- 2 ) 


p = 1.8GeV 


p = 2.6GeV 


q* = 1/a 


g* = 7r/a 


K input 












6.0 


2.44(10) 


123.8(5.7) 


113.6(5.3) 


84.6(3.9) 


80.4(3.7) 


6.2 


1.72(11) 


114.4(3.3) 


110.1(8.0) 


84.5(6.2) 


81.0(5.9) 


6.4 


1.350(57) 


110.3(5.2) 


108.2(5.1) 


86.5(41) 


83.4(3.9) 


a 2 ^0 




104.7(7.4) 


106.0(7.1) 


86.9(5.6) 


84.2(5.4) 


(ft input 












6.0 


2.73(15) 


138.5(8.2) 


127.1(7.5) 


94.6(5.4) 


90.0(5.2) 


6.2 


2.04(23) 


136(16) 


131(15) 


100(12) 


96(11) 


6.4 


1.597(97) 


130.4(8.4) 


128(8.3) 


102.3(6.6) 


98.7(6.3) 


a 2 ^0 




128(11) 


129(12) 


105.5(9.1) 


102.2(8.7) 



TABLE IV. The final results for the strage quark mass m^ /[S (2GeV) obtained with K meson 
or cf) meson mass to fix the bare strange quark masses am s . 
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